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We generalize the surface growth model of Gates and Westcott to arbitrary inclination. The 
exact steady growth velocity is of saddle type with principal curvatures of opposite sign. According 
to Wolf this implies logarithmic height correlations, which we prove by mapping the steady state of 
the surface to world lines of free fermions with chiral boundary conditions. 
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I. INTRODUCTION 

In its simplest form, deterministic surface growth is governed by 

^/i(x,i) =«(V/i(x,t)). (1) 

Here /i(x,t) is the surface profile at time t relative to a reference plane, x S M^, and (|l|) just expresses that the local 
velocity v along h depends only on the local slope u, u = V/i. Usually microscopic growth mechanisms, like deposition 
and evaporation of atoms, are noisy. Kardar, Parisi, and Zhang modelled this randomness by adding noise and 
relaxation to (Q). Expanding in the slope relative to Uo, they obtained the KPZ equation (in 2 + 1 dimensions) 

d , , . -sr^ „ , , (9/i 1 ^ ^ , . dh dh . , ^ /„s 
-/. = «(u„)+^9,z;(u<,)— + - ^ a,:a,„(u„)__ + z.A;i + e (2) 

with ^ space-time white noise. If didjV — Si,jX, A / 0, then the nonlinearity is of the standard form iA(V/i)^. The 
large space-time scaling behavior is then governed by a non-Gaussian strong coupling fixed point. We refer to f^-jj] 
as surveys. 

Wolf |5| noticed that the Taylor coefficient d^v is really a symmetric 2x2 matrix which allows for the possibility 
det d'^v < 0. Based on a one-loop RG calculation he concluded that in this case the nonlinearity in (|^) is irrelevant 
and the statistics of the growing surface is identical to the one of an equilibrium interface. In particular, in the steady 
state the height-height correlations grow logarithmically with distance. 

At the time the prediction of Wolf came as a great surprise. In our paper we support his claim by a model surface 
for which we determine the growth velocity and the steady state correlations explicitly. Our model was proposed 
^ and studied by Gates and Westcott ||] who, in particular, computed the steady state. Since detailed balance is not 
O ■ available, there is no general method for obtaining the steady state. It is defined only implicitly as the invariant 
^ ' distribution of some stochastic evolution. Thus an explicit example is a most welcome exception. Gates and Westcott 
J> , considered only tilts of the form u = {ui, 0), which does not suffice to check the Wolf scenario. To extend their surface 
model to arbitrary tilt and to obtain its steady state is straightforward. One merely has to impose chiral boundary 
. ^ conditions for the height lines. Our main advance is to handle such boundary conditions in the limit of large volume 
' through a transfer matrix which includes the "center of mass" coordinate. We compute the average velocity v(u) 
' ' ' and show that detd^v < for arbitrary u ^ (0,0). To complete the argument we also determine the height-height 
correlations and verify their logarithmic behaviour. Presumably with some extra effort along the lines of Q] one could 
also show that the statistics is Gaussian on a large scale. 
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Let us define the Gates- Westcott model of a crystalline surface (see Figure |^) . At a given time it is specified by the 
height profile h{m, r). Here m G Z, r G M, and h Q Z. If also t £ Z, then the bulk crystal would be built up from 
unit cubes. Mathematically it is somewhat simpler to allow for real r. We assume that h{m, r) is non increasing in 
m. Therefore we can introduce the collection of height lines {^^(t)} such that 

h{m, t) = J - j for 0j (t) < m < 0^+1 (r). (3) 

and (pj+i are the boundaries of the terrace with height J — j. We want to allow for arbitrary step heights along 
the discrete 1-direction. Formally, this means that we admit also a terrace width 0. The height lines satisfy then the 
constraint 0j(t) < 0j+i(r) e Z for all t G M. On the other hand along the continuous 2-direction we only admit step 
height 1. Thus a height line ipj is specified by its location, say at t = 0, (j)j{Q), and the location of the rightward steps 
(=kinks), 4>j{T^) + 1 = (/<j(T_|_), and the leftward steps (=antikinks) , (j)j{T^) — 1 = 0j(T+) as in Figure 




FIG. 1. A piece of the crystalline surface. The deterministic motion is indicated. 

Physically we consider a crystal which is in contact with its vapor phase at a pressure slightly higher than equi- 
librium. Atoms from the gas phase are deposited on the surface. We neglect evaporation. Energetically the fastest 
process is deposition at (anti)kink sites. In first approximation this leads to a constant speed motion of kinks and 
antikinks along the steps. The second fastest process is deposition at steps away from kink sites, which is modelled as 
purely random. Nucleation on the terraces is neglected. For a single height line the dynamics is known as polynuclear 
growth (PNG) model [|. To summarize, the surface dynamics has a stochastic and a deterministic piece. In the 
deterministic part every kink moves with velocity —w, > 0, and every antikink with velocity w along the 2-direction. 
Whenever a kink and antikink pair meet they simply annihilate each other. In addition kink-antikink pairs are ran- 
domly generated with uniform rate ^ provided the constraint 4'j{T) < 4>j+i{T) is not violated. This means, if allowed, 
one has (f>j{To) — > 4'ji'''o) + 1 at some tq. Under the deterministic motion the newly created kink-antikink pair at 
{tq, (/jj^To)) immediately separates with speed 2w. Since the higher lying terraces expand at the cost of lower lying 
ones, the surface grows along the 3-direction. 

In order to determine the steady state we have to work in fin ite v olume with the slope u imposed by appropriate 
boundary conditions. This will be done in Section ||. In Section III we map the height lines to fermionic world lines. 



Their weight corresponds to noninteracing fermions. However there is a complicating constraint due to the M2-tilt of 
the surface. We compute the free energy and the two point function of the constrained fermionic system. In Section 
we relate these results to the growth velocity and the height-height correlations. 
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II. STEADY STATE 



We restrict to periodic surface configurations obeying 4)j+j{T) = 4>j{T) + B and '/'^(t + T) = (/'^(t) + M, where 
T e M+, B,J,M G N. This induces the slope u = (■Ui,it2) with ui = -J/S and U2 = -uiM/T (see Figure 0). 
For convenience we choose ui < 0, M2 > 0. The remaining inchnations foUow by symmetry. In order to determine a 
surface configuration it then suffices to specify J height hnes (/>i(r), . . . , 0j(t) modulo B in the interval < r < T. 
The height in the 3-direction is encoded only modulo J. Let a;["''' G [0, T] be the position of the l-th kink for height line 
j,l = 1,..., 7i(j) + M and j/^'^ G [0, T] be the position of the ^-th antikink, 1 = 1,..., n^^'> . Wc set n = (n^^) , . . . , n^-^)), 

||n|| = 2(n(^) + h n^-^^) + MJ, the total number of kinks and antikinks, and $0 = (0i(O+), . • . ,<l>j{0+)). The 

space of height line configurations F decomposes then as the disjoint union F = Un ^o^^'^' '^O-'' ^^'^'''^ r(n;<i>o) is 
some subset of [0,T]ll"ll as defined through the constraints already explained (we define [0, T]" to be a single point). 
Furthermore for fixed n the F(n; $0) are glued together in a way that is determined by the condition that (anti-)kinks 
leaving the interval [0,T] reappear smoothly on the opposite side, forming the manifold F(n) (we avoid the subtleties 
of defining properly the manifold F(n)). The stationary measure P restricted to F(n) is given by 

^, -Miinii)n n dxpn^^^/'O ' ^(")=^"/^- (4) 
3=1 ^ 1=1 1=1 ^ 



Here 77 = a/^/2w and the normalization 

oc 

^ = E^" E |r(n;'i>o)l (5) 

"=0 n,<E>o;l|n||=n 

with I • I denoting the ||n||-dimensional volume (for ||n|| = we set |F(n;<i>o)| = 1). 

To verify the stationarity of P we first note that the Lebesgue measure on F(n) is left invariant by the deterministic 
part of the dynamics. At a configuration with a total number of ||n|| kinks and antikinks there is a gain in probability 

with rate 2?i;V0p(||n|| + 2) through the annihilation of kink-antikink pairs. Here = lo dT9{(j)j (r) — (/>j-i (r)) , 

6{x) = 1 for a; > and 0{x) = otherwise, is the "number" of configurations immediately before the annihilation event 
which yields (p. The loss in probability occurs with rate ^V0p(||n||) since is also the "number" of configurations 
which arise from by generating a kink-antikink pair. Equating gain and loss we obtain (^. For ?7 ^ the process 
is ergodic on F and therefore the stationary measure P is unique. (Any configuration can be evolved to a reference 
configuration by properly adding kink-antikink pairs in the course of the deterministic evolution). 

The measure P and the partition function Z depend on T, B and J, M. We want to take the infinite volume limit 
B 00, T ^ 00 choosing J, AI such that the slope u is fixed. 



III. FERMI LIQUID WITH CHIRAL BOUNDARY CONDITIONS, FREE ENERGY 

From (^) we see that, subject to the constraint of not crossing, every configuration of height lines has the weight 
7yl|n||_ This is precicely the weight of free fermion world lines in the Euclidian (imaginary time) set-up Q. The surface 
tilt induces chiral boundary conditions. 

Fermion world lines do not overlap. This can be achieved simply by redefining 

^^-(r) =0,(r)+j-l, j = l,...,J, (6) 

with periodized version 

0Jj{t) = Wj(r) mod N, N = B + J , (7) 
(see Figure H). We introduce the relative center of mass 

m{T)= / dr'^d;,(r') G Z. (8) 

Clearly m(-) increases by 1 at every kink and decreases by 1 at every antikink. The chiral boundary conditions are 
then 

m(T) = JM , [bj{T) = [dj{0) + M mod N, j = l,...,J. (9) 



3 




FIG. 2. The mapping of the surface in Figure to fermion world lines. Left and right boundaries are identified. The upper 
and lower boundaries are identified modulo a relative shift of M sites. 

The Fermi liquid lives on {1, . . . , N}. Let {a^, ax} be the corresponding Fermi field operators. They act on the 
fermion Fock space which is spanned by {a+ . . . a+ with jfi) the Fock vacuum. Since m(T) S Z, we introduce 
the auxiliary Hilbert space ^2 with the standard orthonormal basis m G Z}, \m) = (. . . , 0, 1, 0, . . .) at the m-th 

position. On this space let D be the (unitary) right shift, D\m) = \m + 1) . We define the transfer operator 



N 
x=l 



if = -7? V(a+, + a+aa;+i o-n+i 



ai, 



(10) 



acting on !F ^ £2- It expresses that a (anti)kink is linked to the corresponding change in m{T). Let |a;i, . . . , xj; m) 



. |f2) (g) \m) with pairwise distinct xj. Then the matrix element 



.,yj;m,\e '^"\xi, 



(11) 



is given by the sum over all fermion world lines with weight 77", n the number of fermion jumps, and subject to no 
overlap, such that a)j(O) — Xj, (^j{T) = yj, j = 1,..., J, and m(T) = m. For the partition function Z we need 



m(T) = JM and 
defined by 



Then 



M mod N. A fermion configuration is shifted through the total momentum operator 



iP —iP 

e a^e = a^+i 



Z = (0|trj {e'P ®D') 



o-N+i = ai. 



)J\-M^-TH 



|0) 



with tr,/ the trace over the Fock subspace containing J particles. 
By the same argument we obtain for the density correlations 



{{atax){T)a+ay) = Z' 



tr. 



(e'^ D-'ye-^^-^^^'ata^e-^^a+ay] |0). 



(12) 
(13) 

(14) 



In ( p^ and (|^ we take the limit T 00, with M = aT + o(l) G Z, thereby fixing an average drift a G M, and 
subsequently the limit N ^ 00 with J = pN + o(l) G 2N + 1, thus fixing the fermion density p G [0, 1]. We choose 
only odd J to simplify notation somewhat later on. We remark that we expect the final result to be independent of 
the order of limits. However, choices other than the one discussed here look considerably more difficult and wc did 
not persue this issue any further. 

Let us proceed with the partition function. H, P and D are diagonalized through Fourier transform. Let 
A = {l,...,iV}. Then the dual is A* = {-tt + . . . ^ 7^ - 221^ tt}. For Z the dual is [-7r,7r]. We set 
Ofe — N^-'-/^ J2x exp(— ifca;)aa; for /c G A* and represent |m) by exp{—imq) G L^([— tt, tt], (27r)^^dg). Then 



H 



2ri y cos(/c + q)a^ak 



fceA 



(15) 
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considered as multiplication operator in q on L^([— tt, tt], (27r) ^Aq). Also 



fceA* 

on the J-particle subspace. Note that [iJ, e*^ ® D'^] = 0. 
We have 



D'^ = exp - i ^ (fc + q)a^ak (16) 



^ = / § E •'^p (2^7^E + 9) + E(^j + 

'^fei<-</cj ^ 3 = 1 3 = 1 



(17) 



We write cos(fcj + g) = A^^,}. cos(0{fc} + q), where = |e*''i + • • • + e'''''^| and = arg(e'''i + • • • + e*''-^), 

and use the integral representation of the modified Bessel function, /„(x) = tt^^ JJ^ d^/e^'^"'^^'^^'-' cos(n((?!) + q)) for 
arbitrary S M and integer n, to obtain 

Z= GM-iM{ki + --- + kj - J<j){k}))lMj[2T^TA{k}). (18) 

'£l<---<fej 

If ^{fe} is maximal, we call {k} a ground state mode. By periodicity we have exactly N of them. The contribution of 
all other modes in ( [l^ ) is exponentially supressed for large T. Thus 

Z~NlMji2i^TA), (19) 

where A — \ X]k=i exp(27riK/iV)| = iVsin(7r/9)/7r + o(l) for large N. The asymptotic expansion for the modified Bessel 
function, 



In{n/z) = {2mi)-^l^{z + ^l + ^yeMnV'^TT/7^){l + 0{n-^)), (20) 
yields immediately the specific free energy of the fermion system in the infinite volume 

/ = lim lim (TiV)-i InZ = apU/Va + Inry - ln{f], + r/,,)) (21) 



with r/a — anp / {2 simrp) = limT, n-^od M J/2T A and r/^ = ^rj'^ + Va- explained below, rjs ± rja take the role of 
effective right and left jump rates for the fermions. As a further consequence of ( pO| ) we note the relation 

lim = + yiTI^)- exp(-cym^) (22) 

n-^oo In{n/Z) 

for later use. 



IV. CORRELATIONS 

We turn to the expectation values for density, current, and to their correlations. The operator for the density at 
X is px = a^CLx- To obtain the density at imaginary time r we have to evolve a|^aa;with the transfer operator e~'^^ . 
The operator for the rightward fermion current from a; to x + 1 is determined through the limit 

Jx = lim ^(a^+iax+ia^a+)(i)(a+a:rax+ia^+i)(0) 

= V{at+iax ® D) (23) 
by the fermion anticommutation rules. Analoguously the operator for the leftward current from a; to a; — 1, — 

Let us denote ji^^^ = j^, jx''' — VPx and jl — ■ To determine their expectations, we recall that in the limit 
T ^ oo the trace trj reduces to the sum X^ksA* Oki where is the expectation value with respect to the K-ground 
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state a^_7r(/-i)/iv ' ' ' ^t+n{j-i)/N\^)- Because of periodicity these expectation values are invariant under arbitrary 
translations in time, space, and even in dual space. Thus for /3 6 { — 1, 0, 1} we have 

kk' 

= hm Z-'N-'Y [ e2^TAcos(.+9)+zMJ(«+,)g-^(fc+9)/3, + \ 

k K 

= hm iV-iyye-*('=-«)^(a+afe)«Z-i / ^e2^TAcos(K+,)+^(M,/-/5)(K+g) 
T^oo ^ ^ ^ K / / 27r 

= hm N-^Ye-''^^{aXak)^lMj-fiWA)llMjWA) 

T — ^oo ^ — ^ 

k 

= (^aJ/2rjA+ ^lT(^Jj2^y N-'Y,e-'''^{a+ak)o. (24) 

k 

In the limit iV — > oo we obtain 

= P (25) 

and 

UHr)) = i(v/(ap)2 + (2r?sin7rp/^)2 ±ap). (26) 
We define the structure function 

S^^'M = v-'{ji'\r)j^''\0)) - {ji^Hr)){j^''\0)) (27) 
in the infinite volume limit. Assuming r > and using the same method as above, we obtain for finite N 

sf{x,T)^N-^ J2 e-*(^-'^)"-*(^-''+^-''^')((a+a,®i^'')(T)(a+a,. 0i?^')(O)) 

- N-^ e-^('=-'=')^-^('=^+'=''^')((a+afc. ® D'^){T){a+ak ® D'^'){0)) 

kk' 

^-i(k-k')x-i((k+q)l3+(k'+q)fj') 



= hm z-^N-^yy I — 

kk' K 



2r,TAcos(K+g)-2r,r(cos(fe+?)-cos(fc'+9))+iMJ(K+g) / + + \ 

xe [a^af.,a^,ai^)K 



kk' K 

_^g2jjT(A-AA)cos(K-AK+g)+i(Af,/-/3-,3')(K+g) 



27r 



T^oo 

K fcfc 



xe-'^^'-'-P-^'^<^^-hMj-p-p'{2vT{A ~ AA))/lMj{2rjTA) 
= iV-2^e-*(^-^-')--'('=^+'=''^')(a+a,}o(a,,a+}o(aJ/2,7A+ Vl + (a J/2,7A)2)/'+/3' x 



fcfc' 



X exp[-277T^l + (Q;J/277^)2(cosfc - cosk')-iTaJ{smk - sinfc')/A] (28) 

with AA = (T/T)(cos(fc - k) - cosffc' - k)) + 0{N-^T-^), An = {T/TA){sm{k - k) - sin(fc' - k)) + 0{N-^T-^). In 
the last step we used once more (|22| ) and shifted k and k' by k. Taking the limit — > oo results in 

, r^P c\k r'^^^^P Ah' 



-Trp 27r J^rp 27r 



exp[-2T(77s(cosfc - cos A;') + iria{smk - sinfc')) - i(A; - k')x - i{k(3 + k' (3')] (29) 
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for r > and S^^ (x, r) = S^^ {—x, — r) by symmetry. The Fourier transform of S^^ with respect to x equals 

n{K/2)+7rp 

x/ —exp[-ATsiii{K/2){r)sSiiik-ir]aCOsk)-ik{P + P')-iK((3-(3')/2] (30) 

for X > and S^P' {K,t) = S^^^' {-K,t)* for X < 0. For smaU \K\ we have S'^/ {t) oc |X| + 0(iC2), and therefore 
the correlations decay as S^^ (r) ~ for large x and constant r. We remark that the same expectation values and 
two-point correlations are obtained for a pure Fermi liquid with the usual periodic boundary conditions but with a 
non-hermitian transfer operator H = X]a;(^s + '7a)ax+i'^a; + {''Is ~ ila)0'^o.x+i- 

V. SURFACE FLUCTUATIONS 

The growth velocity v of the surface in the 3-direction is proportional to n, the total number of kinks and antikinks 
in{l,...,S}x [0,T], 

V = w{BT)-^{n). (31) 

In the fermion picture n equals the integrated, total leftward and rightward fermion current, i.e. n — 
Eli/o^dr(j+(r) +j--(r)). Thus, using (|2|), 



V = w{B + J)B"V(27?sin7rp/7r)2 + (ap)^ 




because the surface inclination u is related to the fermion density p and drift a as p = |ui|/(l + |mi|) and a = \u2/ui\. 
By symmetry ( |32| ) holds for arbitrarty u G R^. As can be seen from Figure ||, for small u we have qualitatively 
v{u) ~ |u|, which is characteristic for models ignoring nucleation and island formation on terraces. Except for 
u = (0, 0) the growth velocity is saddlelike, i.e. det d^v < 0. Thus the RG analysis of Wolf applies with the prediction 
that the Gates- Westcott crystal surface is in the Edwards- Wilkinson universality class. If so, the surface fluctuations 
should increase logarithmically for large separation. 




FIG. 3. The growth velocity v as a function of the surface tilt u. 
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To check this property for the Gates- Wcstcott model we parametrize the surface relative to the new plane obtained 
by rotating the 1-2-plane by 45° around the 2-axis. The new height function h is then given implicitly for each r by 
h{x, t) = h{m, t) — m where x = h{m, r) + m. h is more explicitly related to the fermion world lines cOx{t), since 

X 

h{x,T) = h{0,T) + Y,^v(.r). (33) 

The height difference between (0,0) and (x,i), A(x,t) — h(x,t) — /i(0,0), equals the number of fermion lines passing 
between (0, 0) and {x,t). We denote the number of fermions at time t in {1, . . . , x} as N[i_x]{t), N[y^x]{t) — YTx'=y Px'{t), 
and the integrated fermion current from site x to site x + 1 during the time interval [0,t] as J[o,t](a;), J[o,(](a;) = 
/q dT(j+(T) — Jx {'''))■ The height difference A{x,t) is then expressed as 

A(x,t) = A^[i,^](0)- J[o,t](a;). (34) 

The mean height difference is given by 

{A{x,t)) ^ p{x - at). (35) 
To compute the height-height correlations we use fermion particle conservation 

J[0,t] {x) - J[0,t] {y) = N[y+l,x] (0) - N^y+l,x] (0 (36) 

for X > y. Then, using that current and density correlations decouple for large distances, we have 

(A(x,i)2)-(A(x,t))2 = 

= ((iV[i,,](0) - J[o,t]ix))^) - (P^ - ^Ptf 

= J™^((^[i,a:](o) - Jmix) + J[o.t]i-y)) X 

x(iV[i,.](0) - J[o,t](2:) + J[o,t](y))) - 
= lim ((7V[-j,+i,,](<) -iV[_^+i,o](0))(%,5;](0) -iV[,+i,^](t))) -pV 

y *oc 

= y'T^ I ^-—^k { . S{k,t)- 



2tt 1 — e^*'^ V 1 — e*'^ 
S{k,0)) 



g-ife(x+l) _ gifc(y+l) ^ s gife(j/— 1) _ g-jfe 



-I 



( ,_,.k sik, 0) ^— , S{k, -t)) 

dk 1 
27r 1 — cos A; 



{S{k,0)-e-'''''S{k,t)), (37) 



with S = 5'°". The rapidly oscillating terms vanish in the limit y oo, since they are multiplied by continuous 
27r-periodic functions. For large keeping x/t constant, the leading contribution is 

(A(x,i)2)-(A(x,i))2 = 

= 7r-2 / — [l-e-^''''^^'''''P cos k{2r]a cos np- x/t) j +0(1) 

= 7r-2ln(|t|) + 0(l). (38) 
If |a;| gets large while t = o{x) we have 

{A{x,t)') - {Aix,t))' = 

= 7r-2y ^ (l-e-^'^l*/"'"' COS A:(2(t/a;)7?a COS 7rp - 1)) +0(1) 

= n-^\n{\x\)+0{l). (39) 

In conclusion we have shown that indeed the height-height correlations increase logarithmically with separation. 
Properly speaking the fluctuations should be considered in the same reference frame as the growth velocity. However, 
this will not change their logarithmic growth. 
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